Exact results of the quantum phase transition for the topological order 
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In this paper a duality between the d — 2 Wen-plaquette model in a transverse field and the d — 1 
Ising model in a transverse field is used to learn the nature of the quantum phase transition (QPT) 
between a spin-polarized phase and a topological ordered state with the string-net condensation. 
The QPT is not induced by spontaneous symmetry breaking and there are no conventional Landau- 
type local order parameters. Instead the string-like non-local order parameters are introduced to 
describe the QPT. In particular, the duality character between the open-string and closed string for 
the QPT are explored. 
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Landau developed a sysmatical theory for classical sta- 
tistical systems. Different orders are characterized by 
different symmetries. The phase transition between dif- 
ferent phases are always accompanied with symmetry 
breaking. However, Landau's theory cannot describe all 
the continuous phase transition such as the transitions 
between quantum ordered phase Topological ordered 
state is a special quantum order with full gapped excita- 
tions. The first nontrivial example of topological order is 
the fractional quantum hall (FQH) states which cannot 
be described by broken symmetries [H- Recently, people 
pay attention to the quantum phase transition (QPT) 
for the topological orders. It is known that there ex- 
ists a quantum phase transition between a spin-polarized 
phase and a topological ordered state. The two phases 
cannot be characterized by a local order parameter. Thus 
such a quantum phase transition cannot be described by 
the Landau theory obviously. To find a new approach 
to learn the nature of such QPT becomes an interest- 
ing issue. The quantum phase transition for the Kitaev 
toric-codc model in a transverse magnetic field was stud- 
ied in Rcf.01. The numerical simulations demonstrate 
the condensation of 'magnetic' excitations and the con- 
finement of 'electric' charges of the phase transition out 
of the topological phase. And then people find that 
the "topological entropy" Stop [1, 0] serves as an or- 
der parameter 0, 0]. At the QPT the "topological en- 
tropy" jumps abruptly from a finite number to zero out 
of the topological phase. In addition, in Ref. Q it is 
revealed that the quantum phase transitions between dif- 
ferent topological orders can be exactly characterized by 
a topological order parameter and shows the one dimen- 
sional (ID) properties. 

In this paper we focus on the Wen-plaquette model 
on a square lattice [l|, Q. The ground states of the 
Wen-plaquette model are Z2 topological ordered state 
which is similar to that for the kitaev toric-code model 
0. After adding the transverse field, a quantum phase 
transition occurs by raising the strength of the trans- 
verse field. However, the model cannot be solved ex- 



actly. Fortunately, we find a duality between the d = 2 
Wen-plaquette model in a transverse field ( transverse 
Wen-plaquette model) and the d = 1 Ising model in a 
transverse field ( transverse Ising model). These results 
indicate a new type of QPT between a topological or- 
dered phase and a non-topological ordered phase, includ- 
ing the non-local order parameters, the ID properties of 
the QPT, in particular, the duality between open-string 
and closed-string. 

The Hamiltonian for the d — 2 transverse Wen- 
plaquette model on square lattice [H, H] is given as: 



(1) 
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and g < 0. af , are 



with Fi 

Pauli matrices on sites, i. 

For the first part of the model, it is an exact solved 
model fg^. The ground state of it is a topological order 
described by Z2 projective symmetry groups. And there 
exists three types of quasiparticles: Z2 charge, Z2 vortex, 
and fermions, of which three types of string operators Tl, 
T2, T3, are defined, respectively The fermions can 
be regarded as bound states of a Z2 charge on even pla- 
quette and a Z2 vortex on odd plaquette. The string 
operator for Z2 charge has a form Wc{C) = Y[ "'''"i 

where C is a string connecting the even-plaquettes of the 
neighboring links, and im are sites on the string. For Z2 
vortex, string operator is Wv{C) = Y[ with C as 

a string connecting the even-plaquettes of the neighbor- 
ing links. For an open Tl (T2) string, the ends points 
represent the Z2 charge ( Z2 vortex). It is known that 
the ground state of the topological order has a conden- 
sation of the closed-strings as (Wc,g(C)) 7^ 0. [l3| On the 
contrary, the open strings have mass gap without con- 
densation. 

It is known that the Wen-plaquette model {h = 0) 
can be mapped onto nearest neighbor Ising chains [11. 
[l2l |. This is because that the Wen-plaquette model ( 
/i = 0) has an energy spectrum which is identical to that 
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FIG. 1: The dotted lines are the dual Ising chains with a as 
the chain-index and i as the site-index. 



of the one dimensional Ising chain llj . Furthermore, we 
find that the transverse Wen-plaquette model on a square 
lattice in Eq.(IT|) corresponds to the one dimensional Ising 
chain in transverse field. In the following part of this 
paper, we will show the duality and use it to learn the 
nature of the QPT between a spin-polarized phase and a 
topological ordered state. 

duality 
commutation 



obtain 
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[Bi, Bj] = 0. Then we denote Ai by t 

B 



and 



we define 
Bi. The two 
relations as: 
[A,Aj] = 0, 
To realize the 



, 

Pauli algebra for Ai and Bj , Bj is denoted by tJ_ i rj_|_ 
Here r^, i , r^, i are Pauli matrices on sites, i + h. Thi 
the mapping between the two models is given as 



i+i ■ 



Then we can explicitly denote the original model by 
the following Hamiltonian, describing the Ising model in 
a transverse field in d = 1: 



a.i~\--. 



(2) 

where a is the chain-index which is different for the d = 2 
transverse Wen-plaquette model on different lattices with 
different boundary conditions (see Fig.l). h is an overall 
energy scale and 2/ = 777 ~ f ^ dimensionless 

coupling constant [l3|. This is a ID Ising model along 
diagonal directions. By the mapping above, we find g- 
term in wen's model corresponds to external field term 
in Ising model, while ft,-term corresponds to Ising-term. 
In addition, an important feature about the duality 



is the boundary condition. The d = 2 transverse Wen- 
plaquette model on N X N square lattice with the open 
boundary condition is dual to 2N — 1 decoupled Ising 
chains ( in a transverse field). Here the chain-index a 
is from 1 to M -|- A*" — 1. However, on N x N square 
lattice with the periodic boundary condition, the d = 
2 transverse Wen-plaquette model 1 corresponds to N 
decoupled Ising chains. 

Global phase diagram — From the duality, one can ob- 
tain the global phase diagram for the original model by 
studying the transverse Ising chain. The model in Eq. ([1]) 
has the same phase diagram to that for the d = 1 trans- 
verse Ising model described by Eq.([2]). Then the original 
model in Eq. ([1]) has a T = quantum phase transition at 
(// = ^ = 1 from topological ordered state with {gj > 1), 
to a gapped spin-polarized state with [gj < 1). For large 
gi, gi > 1, the ground state is topological order which 
is dual to the spin-polarized state in the ID Ising model 
with (t^^^ = . On the other hand, for small gi, 
3/ < 1, the spin-polarized state is in a superposition of 
eigenstates, with different sites uncorrected. This state 
corresponds to the spin ordered state in the dual ID Ising 

model with a spontaneous magnetization ^ 7^ 0. 

The scaling law near the QPT — To obtain the scal- 
ing law near the QPT, the above spin-^ transverse Ising 
model in Eq. ([2|) can be described by the following Hamil- 
tonian 



Hi 



(cL- - Ca,j)(cl,, + i + CaJ+l) (3) 



after 
mation 
fermions IJ, [l 



employing Jordan- Wigner 
of the spin operators 



a,j+i 



transfor- 
to spinless 

- 1, 



2r' r 



a, 3- 



i-iy 



In the fermionic 
Ising chain, the 

- ±2h^ 



lexp(±i7r 4.«Ca,n)(cl ■ + Ca,j). 
n=l ' 

representation for the transverse 
energy spectrum for fermions 

has 



Ek ~ ±2h\J [gi — cos{ka) Y + s\n (ka) has an en- 
ergy gap Ek=o ~ '2h\gj — 1| which tends to zero at 
gi — > 1. So the energy gap for the fermions must tend 
to zero at the transition m = -^^iQi ~ 1) c = 2/ia, 
where a is the lattice spacing. Here m is the mass for 
the Majorana fermion which changes sign by tuning the 
mode across the transition; we have chosen m > to 
correspond to the topological ordered side. [3l • 

The ID properties for the QPT — People have guessed 
that the energy gap for the elementary excitations will 
always close at the QPT for topological orders. Our 
results confirm the conjecture. At the critical point, 
gi = 1, the energy for the fermion excitation becomes 
Ek ~ simc ■ k where the velocity is 2ha. Since the result 
is obtained from the dual Ising chain in Eq.fU), the en- 
ergy gap between the ground state and the first excited 
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FIG. 2: The energy gap 5E (scaled by Ag ) between the ground 
state and the first excited state at the critical point g — h on 
different lattices from analysis results. The inset shows the 
energy gap SE (scaled by <; ) at the critical point g — h on 
different lattices from the exact diagonal numerical results. 
Here N *M denotes a N x M lattice. 



FIG. 3: The relationship between an open-string and a closed- 
string. The hexagon black line denotes a closed-string for 
FiFi+iFi+2Fi+3. The black line below it is a open-string for 
o"f 2) of which two Z2 charges (or Z2 vortexes) live at 
the end denoted by the black squares. 



state is 5Ek ^ c ■ Sk = c ■ where L is the length of 
the Ising chain. The energy gap turns into zero at QPT 
in the thermodynamic limit. In Ref.(0) in the thermo- 
dynamic limit, it is predicted that the gap between the 
ground state and the first excited state closes scales as at 
the critical point SE{J\f) = Ei - ^ M'^ , {M = N"^). 
On &N xN lattice with periodic boundary condition, the 
original model is dual to N decoupled Ising chains, then 
the gap is determined by the length of each Ising chain 
which indeed scales as 5E ^ ~ J\f~^. The Fig. 2 

shows the energy gap 6Ek at the critical point on differ- 
ent lattices from above results and those from the exact 
diagonal numerical results. From the results in Fi^.2, the 
energy gap doesn't always scale as 5E{N) ^ Af~^ . 

In addition, the correlation function at the QPT shows 
the ID characters. The long-distance limit of the static 
correlator for FiFj of the 2D transverse Wen-plaquette 
model can be described by the dual correlator for T^iT^i ji 
of the ID Ising model in the transverse field. For a diag- 
onal CciSG, FiFj — FiFj —ne^ +n ey i thc corrclatioii function 
at the QPT is given as[l^[l6l| 



7r2(4n2 



1 



(4) 



Here to, n are integer numbers. 

The non-local order parameters — Firstly we try to 
use the expectation value for erf and Fi to be the order 
parameters to describe the QPT. However, (erf) and {Fi) 
are finite in both phases but discontinue at the QPT. So 
one cannot use the expectation values of erf and Fi to be 
the order parameters to describe the QPT. 



Instead, we define two kinds of non-local order pa- 

i 

rameters : one is 0i = ( J| erf ); the other is 02 = 

i„ = l 

/ ^ \ 

( n ) ■ Here in = i — nBx + niy represents the sites 
along diagonal direction. Let's explain the physics prop- 

i 

erties of them. On the one hand, Yi ^ Tl or T2 

open-string operator with two ends, z„ = 1 and z„ = i. 

i 

On the other hand, W Fi^ is a Tl or T2 close-string op- 

erator around the points i„ = 1 and in ~ i. Fig. 3 shows 
the relationship between an open-string and a closed- 
string. The expectation values of (open or close) string 

i i 

operators Yi '^f ^^'^ Yi ^"^^ given by the Fig. 4. 

From it one can see that the QPT can be characterized 
in terms of the two kinds of non-local order parameters. 

The open-string - closed-string duality — Another im- 
portant feature for the QPT between the topological 
order and the non-topological order is the open-string 
- closed-string duality. From the duality relationship, 

ffcrL^ af,- erf, - i > t^, i and erf i > i t^, i, 

one can define the first order parameter as (f>i = 



^ 1 I 1 ^ ■ I 1 

a, 1+4 a:«+4 



and the second as 



n 



the dual model. And 01 is just the local order parame- 
ter in dual model which is dual to 02, a 'disorder order 
parameter' for the transverse Ising chain. This prop- 
erty leads to the duality between open-string and closed 
string. 
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FIG. 4: The expectation values of (open or close) string op- 

i i 

erators W af^ and Yl Fin- 



From Fig. 4, one can see that in the topological ordered 
phase, gi > 1, we have {4>i) = and (02) ^ )^ 
This result shows that in this phase the Tl and T2 close- 
strings arc both condensed while the open-strings of them 
arc not. It is consistent to the fact that the ground state 
of the topological order has a condensation of the closed- 
strings {W,,c[C)) ^ 0. 

On the other hand, in the spin-polarizcd phase, gi < 
1, a long range order for dual spin correlation function 
of T^^_i indicates the condensation of the open-string 

operator near the QPT0 



1 ) = lim ( TT erf ) = lim 



It is known that the open-string operator, Yl '^1 



(5) 
will 



generate two Z2 vortexes (or Z2 charges) at both ends 
of the string. The nonzero expectation value of it means 
that it is the 'condensation' of both Z2 vortex and Z2 
charge that breaks down the topological order. The dis- 
appearance of the condensation of closed-string ((^2) = 
or (Wc^siC)) ~ means that this phase is not a topolog- 
ical order any more. 

The situation here is much difference from that in 
Ref.[l|, of which the topological order is broken down 
by the condensation of 'magnetic' excitations (Z2 vortex 
), and the confinement of 'electric' charges (Z2 charge). 
Thus our results indicate there may exist different roads 
to break down the topological order : the condensation 
for both Z2 vortex and Z2 charge or the condensation of 
Z2 vortex ('magnetic' excitation) together with the con- 
finement of Z2 charge ('electric' charge). 

In summary, a duality between the d — 2 Wcn- 



plaquette model in a transverse field ( transverse Wen- 
plaquette model) and the d = 1 Ising model in a trans- 
verse field ( transverse Ising model) are uncovered. By 
the duality, we find that the d = 2 transverse Wcn- 
plaquette model on square lattice undergoes a continu- 
ous phase transition. The quantum phase transition be- 
tween a spin-polarized phase and a topological ordered 
state with the string-net condensation is not induced by 
spontaneous symmetry breaking and there are no conven- 
tional Landau-type local order parameters. Instead, the 
QPT can be characterized in terms of the two kinds of 
non-local order parameters, e.g., the expectation value of 

i i 

( open or close ) string-like operators Y\ erf and W Fi^ . 

i>i = l " i„ = i 

And the topological ordered state are destroyed by the 
condensation of the open-string operators {Z2 charge and 
Z2 vortex). These results indicate a new type of QPT be- 
tween a topological ordered phase and a non-topological 
ordered phase. 
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